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Abstract 

The initial value problem for the Vlasov-Poisson system is by now 
well understood in the case of an isolated system where, by definition, 
the distribution function of the particles as well as the gravitational 
potential vanish at spatial infinity. Here we start with homogeneous 
solutions, which have a spatially constant, non-zero mass density and 
which describe the mass distribution in a Newtonian model of the uni- 
verse. These homogeneous states can be constructed explicitly, and 
we consider deviations from such homogeneous states, which then sat- 
isfy a modified version of the Vlasov-Poisson system. We prove global 
existence and uniqueness of classical solutions to the corresponding 
initial value problem for initial data which represent spatially periodic 
deviations from homogeneous states. 



1 Introduction 



Consider a large ensemble of particles which interact only by the gravita- 
tional potential which they generate cohectively If the particles are de- 
scribed by a time-dependent distribution function f = f(t,x,v) on phase 
space and if we neglect relativistic as well as collisional effects, then the 
time evolution of the ensemble is governed by the Vlasov-Poisson system 

dtf+v-dj-d,u-d,f=o, (1.1) 
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AC/ = 47r/), (1.2) 

where 

p{t,x) := / f{t,x,v)dv (1.3) 

denotes the spatial mass density generated by /, U = U{t,x) denotes the 
gravitational potential, and t>0 denotes time, x E IR^ position, and v G IR^ 
velocity. The same system — with a minus in front of the right hand side of 
Poisson's equation — models a plasma of electrons or ions which interact by 
electrostatic Coulomb forces. Here we are concerned only with the stellar 
dynamic case. 

In the corresponding initial value problem one prescribes the phase space 

o „ 

density at t = 0, i. e. /(0) = / is a given, nonnegative function on IR . How- 
ever, stated in this form, this initial value problem is ill-posed. It has to 
be supplemented with a boundary condition for the potential U at spatial 
infinity. If one wants to model an isolated system, for example a galaxy in 
an otherwise empty universe, one requires that 

Urn U(t,x) = 0, t>0. (1.4) 

x^oo 

o 

This boundary condition — together with some sort of decay of / at infinity — 
is used in all mathematical investigations of the Vlasov-Poisson system and 
related systems that we are aware of, with the exception of systems on 
bounded domains, considered in the plasma physics case. 

In the context of general relativity it is traditional to consider also cos- 
mological solutions. An isolated system can be thought of as a more or 
less localized deviation from flat space, whereas a cosmological solution is 
a deviation from a homogeneous state, usually one with non-zero, spatially 
constant mass density. The present paper investigates this situation in a 
nonrelativistic, Newtonian setting. Our main result is global existence and 
uniqueness of classical solutions for initial data which represent spatially 
periodic deviations from the homogeneous state. For the applications in 
cosmology see [^. Note that very little is known about the existence of 
Newtonian cosmological models without symmetry. The only available re- 
sult is a local in time existence theorem for the Euler-Poisson system [^. 
The present paper proceeds as follows: In the next section we briefly de- 
scribe spatially homogeneous solutions of the Vlasov-Poisson system where 
/ is locally isotropic in v. There are essentially two types of such solutions: 
solutions which exist for all future times and are expanding in a certain 
sense, and solutions which first expand, then contract, and collapse to a 
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singular state in finite time. In general relativity, these would correspond 
to an "open" or "closed" universe respectively. Next we derive the Vlasov- 
Poisson-type system which holds for a deviation from such a homogeneous 
state. In Sect. 3 we fix some notation. A local existence and uniqueness 
result for classical solutions to the initial value problem is established in 
Sect. 4, together with a continuation criterion for such solutions. An im- 
portant ingredient in the proof of global existence is a bound on the kinetic 
energy of the system. Such a bound will be established in Sect. 5. In the last 
section we finally prove that the solutions exist globally in time. In the case 
of an open universe, where the homogeneous background state exists for all 
future times, the meaning of this statement is clear, and, in addition, we 
obtain information on the asymptotic growth of the deviation. In the case 
of a closed universe, where the homogeneous state exists only for a finite 
time, "global" means that the deviation exists as long as the homogeneous 
state. This should not be mistaken for a merely local existence result: it 
might have been thought that the deviation develops a singularity before 
the final "big crunch" of the homogeneous state, but our result shows that 
this is not so. 

To conclude this introduction we give some references to the literature 
where we restrict ourselves to investigations of the existence of global, clas- 
sical solutions. In the case of an isolated system, global existence of clas- 
sical solutions for the Vlasov-Poisson system with general initial data was 
first obtained by Pfaffelmoser His proof was greatly simplified and his 
asymptotic estimates improved by Schaeffer [^] and Horst 0, see also Q for 
a spatially periodic version of these results in the plasma physics case. With 
completely different methods. Lions and Perthame |§ also obtained global 
existence for the Vlasov-Poisson system. A version of this system incorpo- 
rating special relativistic effects was considered by Glassey and Schaeffer 
who showed that solutions of the relativistic Vlasov-Poisson system develop 
singularities in finite time in the stellar dynamic case if the initial data have 
negative energy j^. However, this system has the defect that it is neither 
Galilei nor Lorentz invariant. A fully Lorentz invariant system, consid- 
ered in the plasma physics case, is the relativistic Vlasov-Maxwell system, 
where the question of global existence of classical solutions for general data 
is open. We refer to |11| and the references therein for global results for 
special data. A fully relativistic version with gravitational interaction, the 



Vlasov-Einstein system, was investigated only recently in [12|, and global 
existence was obtained for small, spherically symmetric data. As mentioned 
earlier, all these results treat the case of an isolated system. 
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2 Homogeneous solutions and the system for de- 
viations from these 



In this section we first investigate a class of solutions to the Vlasov-Poisson 
system ( |1.1| ) , ( |1.2[ ) , ( |1.3| ) which are homogeneous in x and locally isotropic in 
V. Then we shall derive a system of equations which is satisfied by deviations 
from such a homogeneous state. 

For a nonnegative function H gC^ (IR) we set 



Ut,x,v):=H(a'{t) 



a(t) 2 
a{t) 



(2.1) 



where a is a positive, scalar function to be determined later. We obtain 



Po{t,x) 

and, after normalizing 



fo{t,x,v)dv = 0-^(1) H{v^)dv 



I' 



H{v^)dv = l, 
we have the homogeneous mass density 

po{t) = a'^{t), t>0. 
A solution of the corresponding Poisson equation is then given by 

Uo{t,x):=—a~^{t)x^, t>0, xGlR^. 
o 



(2.2) 



(2.3) 



This potential does not satisfy the boundary condition (|1.4|) , but this cannot 
be expected from a solution of the Poisson equation ( pT^ ) with a spatially 
constant, non-zero density. We note that 

dxUo{t,x) = —a-'^{t)x, 

and it remains to determine the function a in such a way that /q satisfies 
the Vlasov equation with force term —dxUo- This is the case if the quantity 



a\t) 



a(t) 2 

V —X 

a{t) 
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satisfies this equation. A short computation shows that 
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(dt + v-dx- dxUo ■ dv'j (t) 



ait) 

V rrx 

a{t) 



.2a{t) (v- ■ X (a{t) + —a-\t) 



\ a{t) J V 3 

which is zero for t > 0, x,w € IR"^ if and only if a is a solution of the differential 
equation 

a + —a'^ = {). (2.4) 

O 

This equation describes radial motion in the gravitational field of a point 
mass, and it is well known and obvious that the quantity 

Ea:=\a{tf -^a-\t) 

is conserved along solutions. Let a : [0,Ta[-^]0,oo[ be a solution of Eqn. ( |2.4| ) 
on its right maximal interval of existence with a(0) = 1. We can distinguish 
the following two cases, depending on whether the energy Ea is positive or 
negative: 



Case A, a(0) > a/Stt/S: The solution of Eqn. (2.4) exists on [0,oo[ and 
d{t) > for t>0, i. e. the solution expands forever. 

Case B, d(0) < \/87r/3: The maximal existence time Ta is finite, 
limt^T^ a(t) = 0, and there exists a time ta&[0,Ta[ such that d(t)>0 for 
t E [0,ta[ and a{t) < for t €.]ta,Ta[, i. e. the solution expands up to the time 
ta and then contracts and collapses to a singular state in finite time. 

We remark that in both cases the function a can be extended into the past 
only for a finite time, after which it collapses to a singular state, the "big 
bang" . 

We will now formulate the system of equations which is satisfied by 
perturbations of the homogeneous state fo,po,Uo. Let 

/ = /o + 5, P = Po + ^, U = Uo + W 

be a solution of the Vlasov-Poisson system (|l.l|) , (|l.2| ), (^). A short com- 
putation shows that this is the case if and only if g,a,W solves the system 

dtg + v-dxg- (dxW + ^a-''{t)x) ■d,g = dxW -djo, (2.5) 
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ATy = 47ro-, (2.6) 
a{t,x):= g{t,x,v)dv. (2-7) 



Essentially, one can now consider two types of such deviations: One can 
consider deviations from the homogeneous background which are in some 
sense localized, say, the support of g is compact initially, or one can con- 
sider deviations which extend over the whole universe in some periodic way. 
However, due to the non-zero right hand side in Eqn. (|2.5D an initially com- 
pact support of the perturbation g is not preserved by the time evolution 



of the system (|2.5| ), (2.6), (p. 7]). Also, it seems physically more convincing 
to have the deviation from the homogeneous state extend over the whole 
universe instead of only some part of it. The explicit x-dependence of the 
coefficient in Eqn. ( |2.5| ) seems to preclude the possibility of periodic solutions 
to this system but this can be overcome by the following transformation of 
variables: 

x = a~^{t)x . x = a(t)x 



v — a ^{t)d{t)x ^' ^' v = v + a{t)x^ ^^"^ 



and 



g{t,x,v) = g{t,x,v) = g{t,a{t)x,v + d{t)x), (2.9) 

a{t,x)=a{t,x)=a{t,a{t)x), (2.10) 

W{t,S:) = W{t,x) = W{t,a{t)x). (2.11) 
Then g,a,W satisfy the system 

dtg + -v • - - (d^W + dv) • d,g = 2aH'ia^d^)v- d^W, 
a{t,x):= / g{t,x,v)dv. 



Since throughout the present paper we will be concerned with this system 
in the transformed variables, we now drop the tilde and obtain the following 
system, describing the time evolution of a deviation from the homogeneous, 
locally isotropic state: 

dtg + -v- a,5 - - (d,W + dv) ■ d,g = 2aH'{a\^) vd^W, (2.12) 
a a ^ ^ 

AH^ = 47raV, (2.13) 
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a{t,x):= J g{t,x,v)dv, (2-14) 

where it should be kept in mind that t>0,x,v(z IR^ are related to the original 
variables by the transformation ( |2.8D . We shall be interested in solutions of 
this system which are periodic in x; note that the coefficients in Eqn. ( p.l2| ) 
no longer depend on x explicitly. It should be noted that the resulting solu- 
tions to the original Vlasov-Poisson system are not spatially periodic, since 
we have rescaled the space variable in a time-dependent manner; the change 
in the velocity variable simply means that the bulk velocity a~^{t)a{t)x of 
the homogeneous background solution is shifted into the origin of the veloc- 
ity space. 

To conclude this section, we transform the original Vlasov-Poisson sys- 
tem ( |l . ![) , (|l.2D , (|l.3|) according to (^), drop the tilde and obtain: 

dtf + -v-dxf-- (d^U - ^a-^x + dv) -8,1 = 0, (2.15) 
a a \ 3 J 

A[/ = 47raV, (2.16) 

p{t,x):=J f{t,x,v)dv. (2.17) 

Throughout the paper we make the following assumption: 



, , i a : [O,ra[^]0,oo[ is a maximal solution of Eqn. ( |2.4| ) with a(0) = 1 
^ ' \ HeC^{M), H>0, jH{v^)dv = l, andH{v'^) = 0, |t;|>no>0 

We define fo,Po,Uo by 

27]- 

fo(t,x,v):=H{a^ity), po{t) :=a-'\t), U^{t,x) := —a~\t)x\ (2.18) 



which are the equations (l2.lD , (p. 21) , (^!^), transformed according to (12.81) . 
Then fo,po,Uo is a solution of the transformed Vlasov-Poisson system ( |2.15|) , 
( |2.16| ), (2.17), which we refer to as the homogeneous, locally isotropic back- 



ground state of the Vlasov-Poisson system. If g,a,W satisfy the system 



(2.12), {2J^, KW) then f = fo + g, p = pQ + a, U = Uo + W is a solution of 



t: 



(131) 



2.161) , (|2.17| ), and, via the transformation (|2.8| ), also a solution of 



re original Vlasov-Poisson system ( |1.1| ), ( |1.2| ), ( |1.3D , which we refer to as 
a cosmological solution of the Vlasov-Poisson system, and which no longer 
has to be homogeneous or isotropic. 
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3 Some notation 

Let Q := [0, 1]^ and 5 := Q x IR^. We need to define some spaces of periodic 
functions: 

V{Q) ■= ^h:M^^lR\h{x + a) = h{x), x€lR^, a€2^}, 
p[S) := {h:M^ ^'iR\h{x + a,v) = h{x,v), x,veM^, ae2^}, 

C(Q) := c^{Q)nv{Q), 
:= c"(S)np(5), 

Cl^^^{S) := {heC;^{S)\3u>0:h{x,v)=0, |7;|>u}. 

For [l,oo] we denote by || • \\p the usual L^-norm, where the integral (or 
supremum) extends over Q or 5 as the case may be, and 

LP{Q) :=LPiQ)nV{Q), 
LP{S) :=LP{S)nV{S), 

where we identify functions in V{Q) and 7^(5") with their restrictions to Q 
or S respectively. For the pair we frequently write z. 

4 Local existence and continuation of solutions 

In this section we shall prove the existence of local solutions of the system 
( p^ ), ( pl3D , ( p^ for initial data in Cl^{S). To do this, we have to 
be able to solve the equations ( 2.12| ) and ( 2.13| ) seperately. We start by 
investigating the Poisson equation for a given, spatially periodic density. 

Lemma 4.1 There exists an even function GeC°°{M^\2l^)nV{Q) with 
the following properties: 

(a) Fora^C^{Q) with jQa{x)dx = the function 

W{x):= / G{x-y)a{y)dx= / G{x-y)a{y)dy, x,x'gIR^ 
Jq Jx'+O 



is in C^{Q) the unique solution of 



AW = A-Ka 



with jQW{x)dx = 0. 
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(b) There exists a function Go G C~((IR^ \Z^) U {(0,0,0)}) such that 

G(x) = --^ + Go(x), xe]R3\Z3. 
\x\ 

For a proof of this result we refer to 1^]. Next we note some estimates for 
the potential W which can be proved by treating the singular term in the 
Green's function G as in Props. 1, 2]. 

Lemma 4.2 There exists a constant C > such that for all a as in Lemma 
the corresponding potential W satisfies the following estimates: 

(a) \\W\\oo<c[\\a\ff^\\a\\^J^^ + \\a\\i 

y3|U||2/3. 



or any 



(b) ||a,H^||oo<C(||(7||;/^||^7||^^ + ||a||i 

(c) \\d^,W\\oo<C [{1 + R-'^)\\a\\i + il + lnR/cl)\\a\\oc + d\\dMoc) /< 
d,R withO<d<R<l/2 

(d) \\dlW\\^<C{\\a\\i + {l + \\a\\^){l + \n*\\dM\^/2)), 
where 

s , s < 1 



ln*s: = 



1 + lns , s > 1 



(e) ||a.H^||oo<c(||a||^/^||ar^9 + |k| 



Next we investigate Eqn. ( p.l2 ) for a given potential W . To this end, we 



first consider the corresponding characteristic system 

1 1 

——V, v = -- 

a{t) a{t) 



x = ^v,v = ^—{d^W{t,x) + a{t)v). (4.1) 



Lemma 4.3 Let 0<T<Ta and W eC^'^{[0,T[xQ). Then the following 
holds: 

(a) For z = {x,v) €lR^ and t E [0,T[ there exists a unique solution 
[0,T[3 si-^ {X{s,t,x,v),V{s,t,x,v)) = Z{s,t,z) 
of the characteristic system l \4-ii ) with Z{t,t,z) = z . 
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(b) ZeC\[0,T[^xM^), for every s,t£[0,T[ the mapping Z{s,t,-) isaC^- 
diffeomorphism o/M^, and 

Z{s,t,x + a,v) = Z{s,t,x,v) + (a,0), a € 2^. 

(a) detd,Z{s,t,z)=(^]\ s,tG[0,T[. 
\a{s)J 

Proof : Part (a) and the first two assertions in (b) are standard, the third 
assertion in (b) follows by uniqueness and the periodicity of W . For (c), 
note that 

div^ (-v,--{d^W + av)\=-^-. 
\a a J a 

Thus 

-f-(det9^Z(s,t,z)) = -344deta^Z(s,t,z), 
as ays) 

and the result follows. □ 



Lemma 4.4 Let W satisfy the assumption in Lemma [^.S^ and take 
Cl^{S) such that g{z)+H{v^) > /or 2 G IR^ Then 

g{t,z) := g{Z{0,t,z)) + 2 j\{s)(^H'{a'^v^)v ■d^Wys,Z{s,t,z))ds (4.2) 
= g{Z{0,t,z)) + H{V'^{0,t,z))-H{a'^{ty), te[0,T[, zGlR^(4.3) 

defines the unique solution of Eqn. ^2.1^) in C^([0,r[xlR^) with g{0)=g. 
The solution g has the following properties: 

(a) geC{[0,T[-Cl,{S)), 

(b) Jgg{z)dz = implies Jgg{t,z)dz = 0, t£[0,T[, 

(c) ||5(i)l|oo<||5lloo + 2||i?||oo, t£[0,T[, 

(d) \\g{t)\\i<{\\g\\i + 2)a-\t), tG[0,r[. 

Proof : The fact that Eqn. ( [4.2| ) defines the unique C^-solution of Eqn. 
with 5(0) =g is standard if i?GC2 (IR). Since -H{a^{s)V'^{s,t,z)) is a 



primitive of the integrand in (|4.2| ), Eqn. (|4.3| ) defines the same function, and 
an approximation argument for II shows that it is the unique C^-solution 
also if S C^(IR). The assertion in (a) follows from the periodicity of g and 



10 



Lemma 4.3. To see the assertion in (b), integrate Eqn. (2.12) with respect to 
z, observe that the right hand side in that equation is odd in v, and integrate 
by parts with respect to x and v respectively to obtain the equation 

4: / a{t,z)dz = -'i^^ I g{t,z)dz, te[0,T[; 
dtJs CL{t)Js 

observe that no boundary terms occur due to the compact support in v and 
the periodicity in x respectively. The assertion in (b) follows but observe 
that we obtain no information on 115(^)111 from this since g is not positive. 
The assertion in (c) is obvious by ( |4.3| ). To prove (d), define f-=fo + g- 
Then / satisfies the equation 

dtf + -V ■ d^f --{d^w+dv) ■d,f = o 

a a 

so that 

f{t,z) = f{0,Z{0,t,z)) = H{V\0,t,z)) + g{Z{0,t,z))>0 

by assumption on g, and /(t)GC^(S'), t € [0,T[. The same argument as 
above proves that 



I /(t,z)dz = -3^ / f{t,z)dz, te[o,r[. 

Js a(t)Js 

Thus 



dtJs ' a{t)Js'' 

||/(t)||i = ||/(0)||ia"3(0, 



and since 

||/o(t)!|i = a-3(f), 

the assertion in (d) follows. □ 

We are now ready to prove local existence of classical solutions to the initial 
value problem for the system ( [2T^ ), (gT|), (|2T4|). 

Theorem 4.5 Zei gGCl^^{S) with Jgg{z)dz = and g{x,v)+H{v^)> 
0, {x,v) gIR^. Then there exists a unique maximal solution 

5GCi([0,r[xIR6)nC([0,T[;Ci,,(5)) 

of the system pJ^ j, ^2J^) , with g{0) = g, 0<T<Ta. Moreover, if 

sup{|t;| I (x,?;) G supp5((t), 0<t<T}<oo 

then T = Ta, i- e. the solution exists as long as the homogeneous state. 
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Proof : The method of proof is quite standard, and we only indicate the main 
steps; for an isolated system the corresponding arguments can be found in 

i- 

Consider the following iterative scheme: Define 
and if 

5necH[o,r[x]R6)nC([o,r[;C^,,(S)) 

is already defined, 

an{t,x) := / gn{t,x,v)dv, 



Wn{t,x):=a\t) G{x-y)anit,y)dy 



Q 



for t G [0,Ta[, j;GlR'^, let Zn+i be defined according to Lemma O with Wn 
instead of W, and 

gn+iit,z):=g{Zn+i{0,t,z)) + HiV^^^{0,t,z))-H{a\t)v^),te[0,Tal zGlR^ 

It follows from the above Lemmata that this scheme is well defined. With 

Pn{t) :=snp[\Vn-i{s,0,z)\\z£ supp^U IR^ x (0) , < s < t } 
+uoa-\t), te[0,Ta[, 

it follows that 

gnit,x,v) = 0, \v\>Pnit). 



Using the characteristic system and the estimates from Lemma 4.2 and 
Lemma 4.4 it can be shown by induction that 

Pn{t)<P{t), t(£[0,6[, neM 

where P: [0,5[-^^'^ is the right maximal solution of the equation 

P{t) = Po + C{g) f f P2 (,) + ^P(s) + (s)) ds + u^a'^ (t), (4.4) 
Jo \ a{s) J 

5<Ta. Here C{g) depends only on \\g\\i and H^Hoo, and 

Pq := supjlwl I z G suppr^UlR^ x 
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On the interval [0,6[ the iterative scheme can now be seen to converge. First 
observe that 

\Wnmoo<CPHt), \\d^,W{t)\\^<c[a-\t) + a{t)P\t)), nGiN, te[0,6[. 
Next it can be shown that on any compact subinterval [0,5q] of [0,6[ 

||a.cTn(t)||oo<C, \\dlW{t)\\^<C, neM. 
This can be seen from the following Gronwall-type inequality 

||a.a„(t)|U <Cexp (/^*ln- "^-^--^^"^"- d.] , 

which follows from Lemma (d). Now define 

an{t) := snp{\Xn+i -Xn\{s,t,z) + lK+1 - Vn\{s,t,z) I < S < t, Z G ]R6| 

for tG [0,(^0] ■ Using the uniform bounds which hold on the interval [0,5q] 
another Gronwall argument shows that 

On(i)<C' / an-i{s)ds, 0<t<6o, n>l. 
Jo 

Thus, the sequence Zn is a uniform Cauchy sequence on [0,(5o]^ xlR^. The 
convergence of gn{t), o'nit), Wn{t) in C.,^{S), Ct^{Q), and C^{Q) respectively 
follows, uniformly on [0,(5o]. Using Lemma (c) it can be seen that Wn{t) 
is actually a Cauchy sequence in C^((5), uniformly on [0,(5o]- Thus, the 
function g is , and since < is arbitrary we obtain the desired classical 
solution on the interval [0,(^[. 

To obtain uniqueness, the Gronwall arguments indicated above may be 
used on the difference of two solutions with the same initial datum. 

Let g G C^([0,T[x]R^) be the unique maximal solution, and assume that 
T <Ta and 

P* :=sup{|f I I {x,v) G suppg'(f), < t < T} < 00. 
Then by Lemma |4.4| 

\\9{t)\\i + \\g{t)\\oo<C, te[0,T[. 

Consider now Eqn. ( [4.4D , but with g{to) instead of g and P* instead of Pq. 
Then the function P exists on an interval [to,io + '5[, the length of which is 
bounded from below by a positive constant, which is independent of to & 
[0,r[. Moving to close enough to T the above proof shows that there exists 
a solution on the interval [to,to + (5[ which extends the solution g beyond the 
time T. This is a contradiction, and the proof is complete. □ 
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5 An energy estimate 



Throughout the rest of this paper {g,a,W) denotes a solution of the system 
( |2.12| ), ( p. 13 ), ( 2.14| ) on a maximal interval of existence [0,r[c [0,Ta[ as 



obtained in Theorem 4.5. By C we shall denote constants which may change 
from line to line and may depend on g but never on t. We want to establish 
a bound on the kinetic energy which in turn shall give a bound on ||5(t)||5/3. 
This is desirable because it allows the force term dxW to be bounded by a 
lower power of ||o"(t)||oo, cf. Lemma |4.2| (e). However, the latter argument 
immediately runs into the problem that g is not necessarily positive, and a 
bound on the kinetic energy of g is useless when trying to control ||5(i)||5/3. 
To avoid this problem, we consider f '■= fo + g instead. This function satisfies 
the equation 

dtf + -v-dxf--{dxW + av)-d,f = 0. 
a a 

Thus it is constant along characteristics, and since it is nonnegative initially 
by assumption on g, is nonnegative for all tG [0,T'[. Now define 

Ekinit) ■■= / v'^f{t,z)dz, 
Js 

Epotit) := J^W{t,x)a{t,x)dx = -^-^ J^\dxW{t,x)\^dx, te[0,T[. 

If we differentiate the kinetic energy with respect to time, use the differential 
equation for /, and integrate by parts we obtain 



-Ek^nit) = -5^Ek^nit) + ^ f W {t,x)divj {t,x) dx , (5.1) 

ait) ait) Jo 



d a{t) 2 

^i?..n(t) = -5^£;,.„(t) + ^ 

where 

j{t,x):= / vg{t,x,v)dv; 



/o, being even in v, does not contribute to the mass current density. Inte- 
grating ( p.l2|) with respect to v yields 

1 d 
Ota + -divj + 3-(T = 0, 
a a 

and inserting this into ( ^.1| ) we obtain 

^Ekinit) = -5^ Ekinit) -2 f W{t,x)dta{t,x)dx-6-Epot{t). (5.2) 
at a(tj Jq a 
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Differentiating tlie potential energy yields 
d 



-Epot{t) = 2 f W{t,x)dta{t,x)dx + 2-Epot{t). (5.3) 
at jQ a 

From Eqns. ( f)-2\ ) and (5.3) it follows that 

a^{t)E{t) = E{0)+ (\{s)a^{s)Epot{s)ds, tG[0,r[, (5.4) 
Jo 



where E : = Ekin + Epot denotes the total energy. To exploit Eqn. ( |5.4| ) we 
need to estimate the potential energy in terms of the kinetic energy. This 
argument is due to Horst ^ (5.7)], and we only write down the result for 
our present situation: 

Lemma 5.1 Forte[0,T[, 

\Epot{t)\<c{a~\t) + a-^'/\t)Ek^n{t)'/^). 

To continue, it is convenient to distinguish between Case A and Case B. In 
Case A, d > 0, and since the potential energy is negative, Eqn. ( |5.4| ) implies 
that 

E{t)<a-^{t)E{0). 

With Lemma |5.1| this yields 

Ekinit) < C (a-\t) + a-3/2 {t)Elil (t)) , t g[0,T[, 



which gives a bound on the kinetic energy in Case A. In Case B where d 
eventually becomes negative, we have to argue differently. Observing that 
in this case a is bounded and T < < oo we obtain from (5.4) and Lemma 
5.1 the estimate 

a'^mkinit) < C(^l + J^\d{s)\a\s) [a-\s) + a-^/^Ekinis))'/^)ds 
+ait) + a''/\t)EkUt))'/'^ 
< C^l+ max (a^(s)^fci„(s))^^^^ ^1 + J \d{s)\ds 
<c(l+ max (^a^{s)Ekir,{s) 



1/2 



which yields a bound for the kinetic energy in Case B. Thus we have proved 
the following result: 
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Lemma 5.2 For te[0,T[, Ekin{t) <Ca~''^{t) in Case A, and Ekin{t)< 
Ca~^{t) in Case B. 

By an interpolation argument and estimating g in terms of f + fo, these 
estimates yield the following bounds on ||cr(t)||5/3: 

Corollary 5.3 For te[0,T[, \\a{t)\\^/s <Ca'^^^{t) in Case A, and 
||cr(t) II5/3 <Ca-3(t) in Case B. 

Thus we see that although we have no energy conservation for our problem, 
the usual bounds on the kinetic energy and the consequences thereof still 
hold. 



6 Global existence 

In this section we prove the following global existence result which is the 
main result of the present paper: 

Theorem 6.1 Zei g(^Cl^{S) with Jsg{z)dz = and g{z) + Hiv"^) >0 
for zElR^. Then there exists a unique solution geC^{[0,Ta[xM^)n 
C{[0,Ta[;Cl{S)) of the system ^J^ ), ^^Bj ), ( ^l^j ) with g{0)=g. More- 
over, in Case A, where Ta = 00, we have the estimate 

P(t) = O(t2+'5)^t>0, 

for any 6>0 where 

P{t) :=sup{|t;| I z G supp5(s), < s < t}. 



Proof : The proof makes use of the ideas of Schaeffer 1 13 ] and some unpub- 
lished simplifications of these |14]. Similar estimates have also been used in 
We first consider Case A. Along any characteristic we have the equation 



d_ 

ds 



{a{s)v{s)) = d{s)v{s)-d^W{s,x{s))-h{s)v{s) = -d:,W{s,x{s)). (6.1) 



By Lemma [4.2| (e) and Cor. 



|9x^(5)||oo<Ca2(,)(a-i(s)||a(5)||^9 + ||a(.)||i)<Cia(5)P(.)^/3, (6.2) 
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and thus 

\a{si)v{si)-a{s2)v{s2)\<Cia{t)P{t)^^^\si-S2\, 0<si,S2<t. (6.3) 

Let {X,V) denote a fixed characteristic which hits the support of g, and take 
t G [0,T[ and A € [0,t], T being the length of the maximal existence interval 
of the local solution g corresponding to the initial value g. Then 

\a{t)V{t)-a{t-A)V{t-A)\ 
<[ a^{s)( I [ \g{s,x,v)\dv- 



/x(s)+q7r3 \x-X{s) 

+ ||9a;Go||oo / / \g{s,x,v)\dvdx] 

Jx{s)+QJn^ / 



<C t a-\s)ds + C a\s) f f \g{s,x,v)\dv , 

Jt-A Jt~A JX(s)+QJ]R3 \x-X[S^ 



ds 

dx 
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ds 



where Q:=\— 1/2,1/2]^ . We split the domain of integration of the second 
term into the following parts: 

Ml ■.= {{s,x,v)\D{s,v)<p}, 

M2 := \^{s,x,v)\D{s,v)>pA\x-X{s)\<rD{s,vy'^Y 
M3 := \^{s,x,v)\D{s,v)>pA\x-X{s)\>rD{s,v)'''^Y 

Here D(s,v) ■.= m.m{\v\,\v — V{s)\}, and r>0, p>0 are parameters to be 
chosen later, p<P{t). Let Jj denote the contribution of the set Mi to the 
above integral. The integral Ii is easily estimated: 

* 2/ t - / N dx 



where 



t-A Jx{s)+Q \X-X{S)\^ 



a{s,x) := / \g{s,x,v)\dv 

JD{s,v)<p 

Now 

so that by Lemma [4.2| (e) 



Ii<Cp^/^a{t)A. 
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The estimate I2 is straightforward too: 

2/ \ f f dxdvds 

h< a {s) / ms,x,v) 



2 



t-A Jp<D{s,v)<2P(s)J\x-X{s)\<TD{s,v)-^ ' ' \x — X{s)\ 

< I a^{s) I rD{s,v)~'^dvds 

Jt-A Jp<D(s,v)<2P(s) 

P 

To estimate I3 we want to use the transformation of variables 

z' = Z{t,s,z) 

but since g is not constant along characteristics we switch to f = fo + g first 
and obtain the estimate 

^3 < / / f \f{s,x,v)\lMAs,z)dv- ^-TT^ds 

Jt~A Jx{s)+QJn^ \x-X[s)\'^ 



dx 



+ / a'^is) / „ / \fo{s,x,v)\dv- TT. ..^ ds 

Jt-A Jx(s)+QJn^ \x-X[s)\'^ 



<J + cf a-\s)ds 

Jt-A 



where J denotes the integral over / which remains to be estimated. Let 
A < min{i, 2^^}. Since d is positive and decreasing, 

^a{s2)<a{si)<a{s2), t- A<si<S2<t. 



Using this estimate and Lemma |4.3| in the transformation of variables intro- 
duced above we obtain the following estimate: 



J<C f a\s)j , \f(t,z)\lMAs,Z{s,t,z))- 



dzds 



lzit,s,xis)+Q,n') " ' " ' ''\X{s,t,z)-X{sW 
Now observe that 

Z{t,s,X{s) + Q,1R^) nsupp/(t) C iX{t) + 2Q) x IR^ := M{t) 
if we assume that P(t)A<l/4 i. e. A < min{i, 2^^^, ^p^}. Defining 
J{z):={se[t-A,t]\{s,Z{s,t,z))eMs} 
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we obtain the estimate 

J<C f f{t,z) [ a^{s)\X{s,t,z)-X{s)\-^dsdz. (6.4) 

JM{t) Jj{z) 

To estimate the time integral we need to restrict the lengtli of the integration 
interval so that velocities do not change very much during that interval. We 
have 

a(si) a{si) 

<C,^P{tr/'A + ^Pit)A 
a{si) a{si) 

for si,S2 G [t — A,t]; without loss of generality we may assume P{t) > 1. Thus 
if we choose 

. f 1 1 p 1 

the estimate 

P 

\v{si)-v{s2)\ < - 

holds for si,S2 £ [i — and any characteristic which hits supp^. Now take 
zeM{t) and assume J{z) / 0. Then for se[t — A,t] and a G J{z) the above 
estimate and the definition of the set M3 imply 

^\V{a,t,z)\< \V{s,t,z)\ < ^\Via,t,z)\, (6.5) 

^\Via,t,z)-Via)\<\V{s,t,z)-Vis)\<^\Via,t,z)-Via)\. (6.6) 
Define 

d{s):=X{s,t,z)-X{s), se[t-A,t], 
choose se[t — A,t] such that |(i(s)| = min{|(i(s)| — A,t]} and define 

d{s):=d{s) + d{s){s-s), se[t-A,t]. 

Since d{s) ■ d{s) (s — s) > for se[t — A,t] we see that 

|(i(s)|>|cl(s)||s-s|, s€[t-A,t]. (6.7) 
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On the other hand 



\dis)-dis)\<^\V{s,t,z)-Vis)\ 



^ \d^W{s,X{s,t,z))-d^W{s,X{s))\ 



a2(s) 
and thus 

K.)-<fWI<c.^|»-.1A<l^|.-.| 



By definition of the set M3 and Eqn. (|6.6D 

p<^a{s)m\ 

so that 

\d{s)- d{s) I < ^\d{s) \\s-s\, se[t-A,t]. 
Together with Eqn. (|]^) this yields 

\d{s)\>^\dm\s-s\, se[t-A,t]. (6.8) 



For s G J{z) the definition of the set M3 and the estimates ( |6.5[) and 
imply that 

|d(s)|>Crmin{[t;[>-y(t)|}-l 
Combining this estimate with (|6.8|) and 



\d{s)\ = ^\Vis,t,z)-Vis)\>-^\v-V{t)\ 
a{s) a{s) 

yields 

K.)|>Cmax{-i^|.-y(t)||.-.-|,^,^— ^ 

for sG[t — A,t]. Inserting this estimate in the time integral in ( |6.4| ) and 
distinguishing the cases |f | < |f — and \v\ > \v — V{t)\ finally yields the 
estimate 

a^{s)\X{s,t,z)-Xis)\-^ds < C7a^(t)r->p. 

Jiz) 
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Thus by Lemma 5.2 



J < Cr-^a^{t) / v^f{t,z)dz 

J(X(i)+2Q)x]R3 



< Cr-^a^{t) / v^f{t,z)dz 



Combining the estimates for Ii, I2, and I3 we obtain 
\a(t)V{t)-a{t-A)V{t-A)\ 

< C ( l + a(t)//3 + a2(t)rln^^ +r-iA-i') A 



P 

for 



A = A(t) = min{t,^ 



2a(0)'4P(t)'5C2P(t)4/3. 

Now choose r = a-i(t)P(t)^/2ln"^/2(2P(t)/p) and p = P{tY/^. Then for t 
close enough to T, say t>t*, we have A(t) = 5(g^p(f) and 

\a{t)V{t) - a{t - A)V{t - A) | < Ca{t)P{t)^/^+'A{t) 

for any e>0. Defining tQ:=t and tj+i :=ti — A(tj) as long as ti>t* and 
observing that A(ti) > A(to) we see that there exists A; € IN such that tk+i < 
t*<tk. Thus 

k 

\ait)V{t)\ < aif)\Vit*)\ + Ca{t)P{t)^/^-^'Y.^{U) 

1=0 

< Ca{t){l + t)P{ty/'^+', t(£[t*,T[. 

We conclude that 

P(t)<C7(l + t)P(t)i/2+^ te[0,T[, 
which proves the asymptotic estimate for P{t) and by Thm. shows that 

r=oo. 

Let us now consider Case B: Assume that T <Ta. Then ci<a{t)<C2 
and |a(t)| <C3 for all t€ [0,r[ with positive constants ci, 02,03. Now using 
exactly the same estimates as above we see that 

\a{t)V{t) - a{t - A)V{t - A) | < CP{t)^''^+' A{t) 

for t close enough to T and A(t) : = ^rp^- This yields the same estimate for 
P{t) on the interval [0,r[ as in Case A, in contradiction to Thm. 4.5. Thus, 
T = Ta also in Case B, and the proof is complete. □ 
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